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Abstract 

The choice of a factorization scheme suitable for Monte Carlo simula- 
tions of NLO initial state parton showers is discussed in this contribu- 
tion. 



1 Introduction 

Generating initial state parton showers (IPS) in hadronic collisions at the NLO accuracy is a task 
for which no satisfactory solution has so far been found. An attempt at solving this problem 
which is presented in this contribution is based on an exploitation of the freedom in the choice of 
the factorization scheme. 

As a simple illustration, consider a non-singlet nucleon structure function Fns {x, Q^) . Its 
Mellin moments are given as the product 

i^NsKQ') =CNs(n,^,Fs)gNs(n,M,FS), W 

where CnsC"-; Q/^: FS) stands for Mellin moments of the corresponding coefficient function 
and q^Q,{n,M,YS) represents Mellin moments of the relevant non-singlet parton distribution 
function. Both C^'g,{n, Q/M, FS) and q^si'n, M, FS) depend on a particular factorization scheme 
FS and on a factorization scale M, however, their product ([T]) is independent of both of them. 
The coefficient function Cns(?^, Q/-^, FS) is calculable within the framework of perturbative 
QCD and can thus be expanded in powers of the QCD couplant a = Os/vr 

CNs(n,^,Fs) =f;a^(/U,RS)4';](^n,-^,FS,;U,Rs) . (2) 

Although both the couplant a(/i, RS) and the coefficients C^^{n, Q/M, FS, /U, RS) depend on 
a particular renormalization scheme RS and on a renormalization scale /i, which is in principle 
different from M, the series, if summed to all orders, is independent of both the RS and /i. The 
non-singlet parton distribution function gNsC'^j Af, FS) satisfies the evolution equatioij^ 

^^^^^^^^^ = a(M)PNs(n,M,FS)(?Ns(n,M,FS), (3) 



From the relations js) and ji}, we see that the non-singlet distribution function qisis(n, AI, FS) also depends on 
the renormalization scheme in which the renormalized couplant a[M) is defined. This scheme can in principle be 
different from that used for the expansion of the coefficient function C-^sin, Q/M, FS), but usually these schemes 
are chosen to be identical. 



where the non-singlet sphtting function Pns(^; FS) can be expanded in powers of a(M) 



PNs(n, M, FS) = ^\^) 4s FS) . (4) 

fe=0 



In the next-to-leading order approximation, we retain only the first two terms in the ex- 
pansions Q and Q: 



CNs(n,-|,Fs] = cig(n) + a(/i)C«('n,-|,FS), (5) 



PNs(n, M, FS) = P^^l (n) + a{M) P^^l (n, FS) (6) 



NS V'V ^ ) ^ NS 

and the NLO couplant a(/u) obeys the differential equation 

da(/i) 



dln/i 



6a2(/x)(l + ca(/i)). (7) 



The LO terms CN'g(n) and pj^^{n) are universal — independent of any unphysical quanti- 
ties such as renormalization and factorization scales and schemes. The NLO contributions 

Cjjg^(n, Q/M, FS) and Pj^^iin, FS) satisfy the following condition 



4sYn,|,Fs')=<)(n) 



K(n)+P(5(n) In ^ + ip(^s)(n,FS) 



(8) 



where ^(n) is a scale and scheme factorization invariant. The ambiguity related to the fac- 
torization procedure is already at the NLO large, because the splitting function Pj^g'(x,FS) is 
completely arbitrary function — for any function /(x), there always exists such a factorization 
scheme FS in which P^^{x,FS) = f{x). There are two prominent choices of the splitting 
function P^^{n, FS), which are in some sense opposite to each other. In the first one, the split- 
ting function Pj^g''(n, FS) is set equal to zero. For this choice, which will be called the ZERO 
factorization scheme, the evolution of the distribution function gNs(^) M, FS) is formally iden- 
tical to the LO one and all NLO corrections are thus contained in the NLO coefficient function 
C^^{n,Q/M,FS). The latter choice consists in selecting the splitting function Pj^g''(n, FS) in 

such a way that the NLO coefficient function C^^{n, Q/M, FS) vanishes for M = Q (see the 
equation (jsjl). In this case, the relation between the structure function F^s(^x, Q^) and the dis- 
tribution function qt^sin, Af, FS) has the same form as in the LO and all NLO corrections are 
included in the evolution of the distribution function q^f^sin, Af, FS). This type of choice is the 
essence of the so called DIS factorization scheme introduced in [1 ], which is widely used in phe- 
nomenology. The factorization scheme dependence of NLO theoretical predictions for physical 
quantities is studied, for instance, in [2], where only factorization schemes interpolating between 
the DIS and MS factorization schemes are considered. 

The factorization scheme specifies the way in which the so called coUinear singularities, 
which are contained in cross-sections at parton level, are absorbed into the dressed parton distri- 
bution functions. Within the framework of dimensional regularization, the relation between the 



dressed and bare distribution functions is given in the general case by the formula 



+ a^(M)(^A;f(y) + l4^)(,) + Alf(,)) + 



(9) 



The matrices A^^^^\x) can be chosen arbitrarily and their choice fully specifies the factoriza- 
tion scheme. The factorization scheme can also be specified by higher order splitting functions 
p//^(x,FS), A; > 1, which we can choose at will. The most widely used factorization scheme is 
the so called MS factorization scheme, which is defined by setting the matrices A^^^\x) equal 
to zerc|3and is thus convenient for theoretical calculations. 

At present time many QCD cross-sections at parton level are known at the NLO accuracy. 
However, the algorithms that are used for their incorporation in Monte Carlo event generators 
attach to them the IPS only at the LO accuracy because generating IPS at the NLO accuracy is 
very complicated in the standard MS factorization scheme. The reasons for that are basically 
two: The NLO splitting functions no longer correspond to basic QCD vertices and the splitting 
functions at the NLO approximation are negative for some x, which prevents us from using 
straightforward probabilistic interpretation. Because the IPS induce the scale dependence of 
parton distribution functions, it is inconsistent to attach the LO IPS to NLO QCD cross-sections, 
which include NLO parton distribution functions. This deficiency could be removed by choosing 
the ZERO factorization scheme, in which the NLO IPS are formally identical to the LO ones 
and all NLO corrections are thus put into hard scattering cross-sections. The main advantage 
of this approach is the fact that the existing algorithms for parton showering and for attaching 
parton showers to NLO cross-sections need not be changed. The only step necessary to do is 
transforming parton level cross-sections from the standard MS factorization scheme to the ZERO 
factorization scheme and determining parton distribution functions in the ZERO scheme. 



2 The transformation of hard scattering cross-sections 

In the case of a hadron collision, a cross-section (t{P) (in general differential) depending on 
observables P is given by the formula 

^(^) = 5^/ / dxidx2 AM(a;i,M,FS)Z)j/B(x2,M,FS)cTij(xi,X2;P;M,FS), (10) 
Jo Jo 

where Di/y^{x, M, FS) and A/g(x, M, FS) are the parton distribution functions of the colliding 
hadrons. The partonic cross-section aij{xi, X2; P; M, FS) can be expanded in powers of the 
QCD couplant a(^): 

aij {xi,X2; P; M, FS) = ajf (xi, X2; P) + a(/u) alf{x^,X2; P; M,FS) + 0{a\fi)). (11) 

^with the renormalized couplant a{M) defined in the MS renormahzation scheme 



The LO cross-section (tIj\xi, X2; P) is independent of the factorization scale and scheme. 

The dependence of the NLO cross-section alj\xi,X2; P; M, FS) on the factorization scale and 
scheme is determined by the formula 



alj\xi,X2;P;M,FS) = a'ij\xi,X2;P;Mo,MS) + Yl / ^2/ 

k 



ail\xi,yx2;P)x (12) 



What we need for the conversion to the factorization scheme FS is the knowledge of the corre- 
sponding matrix function Tj-j\x, FS) |^ which is process independent. In the space of Mellin 
moments, the matrix Tj^j\n, FS) satisfies the following matrix equation: 

T(^)(n,FS),p(°)(n)j - 6T(^)(n, FS) = p(^)(n, MS) - p(^)(n, FS) . (13) 
3 The ZERO factorization scheme 



The solution of the preceding equation ( 13 1 for the ZERO factorization scheme reads 



Tallin) = T^Un) = T.^\n) - (n,MS), (n) = (n) = Ti'\n), 

t!uIM) = T^IM) = ^i'^N - -AlM^^)^ Tallin) = r« (n) = 4'\n), (14) 
4gN = -^4gKMS) - 2niTP{n), 

where the unknown functions t/^'* (n) can be expressed in terms of the Mellin moments of the LO 
and NLO splitting functions. The Mellin inversion of Tj^^\n) has to be calculated numerically. 
This was performed for three and four effectively massless flavours. The obtained results are 
very surprising because for x <0.1 

Tj;^\x) PS C,x~^ with ^(nf = 3) = 4.63 and ^(nf = 4) = 3.85 (15) 

and the coefficients d are so large that the functions Tj^^\x) strongly dominate over the NLO 
splitting functions P^;^^(x,MS) in this region (x < 0.1). So the question of applicability of the 
ZERO factorization scheme arises. 

The parton distribution functions in the ZERO factorization scheme are plotted for x € 
(10^'^, 10^^) in Figure [1] In this region they behave like x^^ with ^ close to that in the rela- 



tion ( 15 I. The ZERO factorization scheme can thus provide reasonable predictions only if large 
cancellation between positive and negative values occurs in expressions for physical quantities. 
Hence, within the framework of numerical computations, the ZERO factorization scheme can 
only be used in kinematic regions where x > 0.1. The only exception is its application in the 
non-singlet case, where the functions T^*^^^ (x) do not appear and no problems with large numbers 
occur. 



^7;^^' {x, FS) = - ^i"'' {x, FS), where the matrix {x, FS) is defined in the relation 



Quark singlet distribution 



Gluon distribution 




Fig. 1: Comparison of tlie ZERO and MS parton distributions at M = 50 GeV. Tlie MS distributions were obtained 
by evolving the starting distributions of the MRST98 set (Sj with the fixed number of active flavours rif = 3 (only 
light flavours are taken into account). The ZERO distributions were calculated from the MS ones by using a numerical 
transformation method based on Mellin moments. Note that the gluon distributions are plotted in their absolute value 
because the ZERO gluon distribution is negative in the displayed region. The zero point where the ZERO gluon 
distribution changes the sign is close to a; = 0.1. 



4 Summary and conclusion 

The ZERO factorization scheme is optimal for Monte Carlo simulations of NLO initial state 
parton showers. However, because of the problems with large numbers, this scheme has too 
little range of applicability in numerical calculations. The ZERO factorization scheme should 
thus be replaced by some "almost ZERO" factorization scheme which is sufficiently close to the 
ZERO factorization scheme and is free of problems with large numbers. Searching for such a 
factorization scheme has already been started. 
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